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Abstract 

cn ; 

The confluent algorithm, a degenerate case of the second order supersymmetric 
quantum mechanics, is studied. It is shown that the transformation function must 
asymptotically vanish to induce non-singular final potentials. The technique can be 
used to create a single level above the initial ground state energy. The method is 
applied to the free particle, one-soliton well and harmonic oscillator. 

c 

1 Introduction 

> ' 

The second order supersymmetric quantum mechanics (2-SUSY QM), which involves 
second order differential intertwining operators [TJ |21 El EH El EH HI IE] nas proved useful to 
surpass the difficulty of 'modifying' the excited state levels inherent to the standard first 
order 1-SUSY QM. In fact, in the 2-SUSY treatment we do not respect the restrictions 
imposed by 1-SUSY on the transformation functions Ux(x), U2(x), i.e., they can have nodes 
but induce a non-singular 2-SUSY transformation OEI1EJE]- In this way, potentials 
with two extra bound states above the ground state energy of the initial Hamiltonian H 
(or above the lowest band edge if V(x) is periodic) have been recently generated Plll|ll2| . 
A different atypical method employing two complex conjugate factorization energies has 
been as well implemented, generating in this case families of real isospectral 2-SUSY 
partner potentials of V(x) [T3] . 

There is yet another situation worth studying in detail, namely, when the two factorization 
energies tend to a common real value e. This so-called confluent algorithm was used to 
generate a particular family of isospectral oscillator potentials [13] . However, we have not 
detected a generic analysis (for arbitrary factorization energies e) of the properties of the 
transformation function u(x) ensuring that the final potential will be non-singular. This 
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is the subject of the present paper, which has been organized as follows. In section 2 an 
alternative view of the confluent 2-SUSY algorithm will be elaborated. The restrictions 
imposed onto u(x) in order to obtain non-singular final potentials will be analysed in 
section 3. In section 4 we will apply the technique to the free particle, one-soliton well 
and standard harmonic oscillator. 



2 Second order supersymmetric quantum mechanics 

The second order supersymmetric quantum mechanics (2-SUSY QM) is a particular re- 
alization of the standard supersymmetry algebra with two generators 



{Qj, Qk} = S jk H 8S , [H ss , Qj] = 0, j, k = 1,2, (1) 
where Q x = (Qt + Q)/y/2, Q 2 = (Q f - Q)/{i\/2), 



H _ ( AA] ^ _ f (H - e^H - e 2 ) \ ( ) 

ss V A1A J V (H- ei )(H-e 2 )J' [ ' 

and H, H are two intertwined Schrodinger Hamiltonians: 

HA = AH, (4) 
H=-& + V{x), H = -£, + V(x), (5) 
A =£+V(x)£+J(x). (6) 

The relations between r)(x), j(x), V(x) and V(x), are: 

V = V + 2 V ', (7) 
7 = d - V + r] 2 /2 - rf/2 , (8) 
W - iv') 2 / 2 + ri 2 (v 2 /2 ~ 277' - 2F + 2d) + 2c = , (9) 

where, in terms of c € 1, del, the factorization energies read e\ = d + y/c, e 2 = d — \/c. 
Suppose that V(x) is a given exactly solvable potential; then V(x) will be effectively 
determined if we find explicit solutions t](x) to the non-linear second order differential 
equation @. Depending on the sign of c, two essentially different cases arise. 

If c then e\ ^ e 2 , and we must look for solutions of the two Riccati equations: 

P[ + P 2 = V-e t , z = l,2. (10) 

Having 0i(x), ^(x), we get two different equations for rj(x) (see, e.g., jEHUl) 

i = f] 2 + 2/3 lV + e 2 - e x , (11) 
r/ = ?? 2 + 2/5 2 ^ + ei - e 2 . (12) 
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By subtracting them, we arrive to a finite difference algorithm for 7](x): 

V(x) = (e l -e 2 )/(i3 1 -i3 2 ). 



(13) 



On the other hand, the confluent case arises for c = implying that e\ = e% = e = d. In 
this situation, we look for solutions to just one Riccati equation [TJ] 

p , + p 2 = V-e, (14) 

and rj(x) satisfies an equation arising when e\ = €2 in (jllll2j) : 

7/ = + . (15) 

This is the Bernoulli equation, whose general solution is given by: 

r)(x) = —w'(x)/w(x), (16) 

where 

w {x) = w - J e 2 f^ x)dx dx . (17) 

In the next section we will analyse the conditions which grant that the confluent 2- 
SUSY transformations are non-singular. This provides the simplest way of ensuring that, 
departing from an exactly solvable initial potentials V(x), we arrive as well at an exactly 
solvable regular V(x) (see e.g. the discussion in [T5|). 



3 Confluent non-singular transformations 

Let us express first the confluent formulae of section 2 in terms of solutions of the initial 
Schrodinger equation obtained from |T4|) by the change (3{x) = u'{x)/u{x): 

-u"(x) + V{x)u(x) = eu{x). (18) 

Thus, up to an unimportant constant factor (see (jlfij) ). the key function w(x) becomes 

w(x) =w - u 2 {y)dy, (19) 

Jxo 

and the confluent 2-SUSY potential V(x) is given by: 

V(x) = V(x) - 2[w\x)/w(x)}'. (20) 

It is clear now that in order to arrive at real non-singular potentials V(x) we have to use 
real solutions u(x) of (|TS|) inducing a nodeless w(x). Let us notice that 

w \x) = -u 2 (x), (21) 
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meaning that w(x) is decreasing monotonic, so the simplest way of avoiding its zeros is 
to look for the appropriate asymptotic behaviour for u(x). Two different situations are 
worth considering. 

(i) Suppose first that e = E m is one of the discrete eigenvalues of H and the transformation 
function is the corresponding normalized physical eigenfunction, u(x) = ip m (x). Denote 
by v + the following finite integral: 

/>oo 

v+= u 2 (y)dy. (22) 

Jxo 

It is straightforward to show that: 

lim w(x) = wo — u + + l = v + l (23) 

X — > — oo 

and 

lim w(x) = v. (24) 

x^+oo 

It turns out that w(x) is nodeless if either both limits are positive or both negative, 
leading to the zz-domain where the confluent 2-SUSY transformation is non-singular: 

i/eR\(-l,0) = (-oo,-l] U[0,oo). (25) 



(ii) Suppose now that the transformation function u(x) is a non-normalizable solution of 
(HHJ associated to a real factorization energy e ^ Sp(if) such that 



POO 

lim u{x) = and u + = / u 2 (y)dy < oo. 



(26) 



If this is the case we can show that: 

"X 



lim w(x) = wo + I u 2 (y)dy = oo (27) 

x^— oo / „ 

J — oo 

and 

lim w(x) = Wq — v + = v. (28) 

x^+oo 

By comparing both limits and taking into account that w(x) is decreasing monotonic, it 
turns out that w(x) is nodeless if 

v > 0. (29) 
Let us notice that the same //-restriction holds in the case when 

lim u(x) = and is = j u 2 (y)dy < oo, (30) 

oo 



X— > — oo 



though now v = —(wo + z/_). 

Once the regularity of the confluent 2-SUSY algorithm is assured, we analyse the spectrum 
of H. From the intertwining relationship (jlj) and the factorizations in (JHJ) we immediately 
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obtain normalized eigenstates \ip n ) of H provided that v satisfies either ()25|) in case (i) or 
(JUJ) in case (ii) and A\ip n ) ^ 0: 

$ n ) = (E n -e)- 1 A\1> n ), (31) 

(so in the case (i) we cannot obtain \^ m ) of (pITj) because A\ip m ) = 0). The orthonormal set 
{\ip n ),n = 0, 1, 2, . . . } so constructed is not automatically complete (we have to analyse 
yet the existence or not of an extra normalizable function ip belonging to the Kernel of A' 
which is orthogonal to all the \ip n ),n = 0, 1, 2, . . . ). To find ip explicitly, let us factorize 
as follows: 

A=[£ + Pi*)][£-(K*)-v(*)]- (32) 

It turns out that the ip E Ker(A^) we are looking for is annihilated by the second factor 
operator of (pT2"j) : 

${x) = n Q e$ [ ^ x)+ ^ x)]dx = n u(x)/w(x), (33) 

no is a constant. It is straightforward to check that if)(x) is a normalized eigenfunction 
of H with eigenvalue e for v E R \ [—1, 0] in the case (i) with n = a/ v{y + 1) and for 
v > in the case (ii) with riQ = y/u. On the other hand, ip(x) becomes non- normalizable 
for v = —1,0 in the case (i) or for v = in the case (ii). Thus, when e = E m and 
u(x) = ip m (x) it turns out that Sp(if) = Sp(if) if v G R\ [—1, 0] while the level E m is not 
present in Sp(if) for v = —1,0 (in this been 'deleted' in order to generate 

H). On the other hand, when e Sp(if) and u(x) obeys either (j2EJ) or (|3T)|) it turns out 
that Sp(-fiT) = {e} U Sp(H) for v > but Sp(#) = Sp(F) for ^ = 0. For all the other 
^-values [y E (—1, 0) in case (i) and v < in case (ii)) it arises a singularity in V(x) due 
to the existence of a zero in w(x). We notice, in particular, that the case (ii) allows to 
generate a single level above the ground state energy of H, a mechanism which cannot be 
directly implemented in the 1-SUSY treatment. 

Let us remark that our confluent 2-SUSY procedure coincides with the Abraham-Moses 
generation technique of creating, deleting or changing the normalization of a single en- 
ergy level PHI (see also |T71 ll8|). The same procedure, known as binary Darboux trans- 
formations [T§|, has been employed to generate bound states embedded in the contin- 
uum [Hlinil221- 



4 The simplest applications 

Let us analyse now the simplest applications of the confluent 2-SUSY algorithm. 

(a) Consider first the free particle for which V(x) = 0. For a fixed arbitrary e < 
which does not belong to Sp(if ) there are two asymptotically vanishing transformation 
functions: 

u{x) = V2ke ±kx , e = -k 2 . (34) 
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A direct calculation leads to: 

w (x) = T 2e ±k(x+Xl) cosh[A;(x - Xl )}, (35) 

where v = e ±2kxi > 0. By substituting those expressions into (|20[) we obtain the Poschl- 
Teller potential in both cases: 

V(x) = -2k 2 sech 2 {k(x - an)], (36) 
which has a bound state at e = —k 2 . 

(b) Take now the previous Poschl- Teller as the initial potential: 

V(x) = -2k%sech 2 (k x), (37) 
and denote the ground state energy as usual, E = —k 2 . 

Let us consider first the case when e = E and u(x) is the normalized ground state: 

u (x) = y^sech(fc x). (38) 
A straightforward calculation leads to: 

w(x) = v + g — \ tanh(fco#), (39) 
which produces once again the Poschl- Teller potential: 

V(x) = -2klsech 2 k {x - x x ), (40) 
where tanh(/c ^i) = 1/(1 + 2^)- 

Suppose now that e = — k 2 ^ Eq, k e R. The solutions with the right asymptotic 
behaviour are here: 

u(x) = V2ke ±kx [k tanh(k x) T k], (41) 

leading to: 

w (x) = T{v + e ±2kx [k 2 + k 2 T 2kk tanh(k x)]}. (42) 
It turns out that the confluent 2-SUSY potential V(x) acquires the Bargmann form: 

~ , , _ 2{k\ - k\)\k\^&?k x {x - x x ) + /c|csch 2 fc 2 (a; - x 2 )\ , . 

\k\ tanh k±(x — x\) — hi coth k 2 (x — x 2 )] 2 

where for k > k we need to take k\ = k , k 2 = k, v = (k 2 — kl)e ±2kX2 , e ±2k ° Xl = 
(k + ko)/(k — ko) while for k < k we require k\ — k, k 2 = k , v — (k^ — k 2 )e ±2kxi , 
e ±2koX2 = (k + k )/(k - k). 

(c) Finally, let us analyse the harmonic oscillator potential: 

V(x) = x 2 , (44) 
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which has a purely discrete spectrum composed of E n = 2n + 1, n = 0, 1, . . . and eigen- 
functions given by: 

Mx) = {V^Tn^ l ' 2 e- x2 l 2 H n {x), n = 0, 1, . . . (45) 
where H n (x) are the Hermite polynomials. 

Let us suppose first that e = E m with m fixed, u(x) being the corresponding normalized 
eigenfunction if) m (x) of (j43|) . The calculation of (|T5|) with x = leads to: 

mo 

w (x) = „ + I - xE J^^g^ggg^ 2 F 2 ( mi , mi _ s; 5 + I, mi + l - s; -x 2 )(46) 

where 2-^2( a i, °2! &i, z) is a generalized Hypergeometric function [231, m o — ( m — 
<5)/2, mi = (m + 8 + l)/2, 5 = if m is even but 5 = 1 if m is odd. It turns out 
that V"(x) is isospectral to the oscillator potential, a case illustrated in figure 1 for e = 7 
(m = 3) and v = —5/4. Let us notice the already involved explicit expression of w(x) in 




-4 -2 2 4 



Figure 1: The confluent 2-SUSY partner potential (black curve) isospectral to the oscillator (gray curve) 
generated by employing the normalized eigenfunction (|45|l for n = 3, e = £3 = 7 and ^ = —5/4. 

In turn, for e ^ Sp(if) the asymptotically vanishing Schrodinger solutions become: 

u(x) = e-T [^(1=5, |; x 2 ) ± 2xggj ^(2=2, |; x 2 )] , (47) 

where \Fi(a, c; z) is the Kummer hypergeometric series. The explicit expression for w(x) 
is too involved to be shown here (three infinite sums of kind (|4fi|) arise in this case). 
Alternatively, we performed a numeric calculation of V(x) for e = 8 taking the solution 
u(x) of (|4~7j) with the upper plus sign and = —5, x = in (JT9j) (see figure 2). The 
spectrum of V(x) is composed of the oscillator eigenenergies E n = 2n + 1, n = 0, 1, . . . 
plus a new level at e = 8. This illustrates clearly the possibility offered by the confluent 
2-SUSY algorithm of creating one single level above the ground state energy of H. 
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Figure 2: The confluent 2-SUSY partner potential V(x) (black curve) of the oscillator (gray curve) 
generated by employing the Schodinger solution (|47|) with the upper + sign and e = 8, Wq — —5, x a = 0. 
The potential V(x) has an extra bound state at e = 8 compared with the oscillator spectrum. 

We conclude that the second order supersymmetric quantum mechanics is a powerful tool 
for designing in a simple way potentials with given spectra, a subject supplying us of 
solvable models with possible applications in the physical sciences (see e.g. |2*4"]). 
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